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A computer  code,  RAD0T,  Has  been  developed  to  track  Incident  radar 
through  forests  consisting  of  realistically  simulated  trees. 

RAD0T  fires  parallel  radar  rays  at  the  simulated  forest  and  records 
Information  for  all  materials  Intersected  between  the  origin  and  the 
ground.  This  Information  Includes  type  of  object  (branch,  twig  or 
leaf),  Its  size  and  orientation  parameters  (components  of  normal  or  - 
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^axis  vector)  and  the  depth. 

Radot  has  been  successfully  Installed  on  the  ETL  CDC-6L00  computer.' 
Other  versions,  for  IRM-36O  or  CDC-66OO  computers  are  also  available. 
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A computer  code,  RADOT,  has  been  developed  to  track  Incident  radar 
through  forests  consisting  of  realistically  simulated  trees. 

RADOT  fires  parallel  radar  rays  at  the  simulated  forest  and  records 
Information  for  all  materials  Intersected  between  the  origin  and  the  ground. 
This  Information  Includes  type  of  object  (branch,  twig  or  leaf),  Its  size 
and  orientation  parameters  (components  of  normal  or  axis  vector)  and  the 
depth. 

RADOT  has  been  successfully  Installed  on  the  ETL  CDC-6400  computer. 
Other  versions,  for  IBM-360  or  CDC-6600  computers  are  also  available. 
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I . INTRODUCTION 


The  distinctive  radar  "signature"  of  a military  object  Is  a prime  mech- 
anism for  Its  acquisition  as  a target.  A major  problem  tn  target  recognition 
Is  the  radar  background  generated  by  a masking  environment.  This  could  be 
artificial,  such  as  camouflage  nets,  or  natural,  such  as  vegetation,  tn  order 
to  cope  with  this  problem,  It  14  Important  to  be  able  to  quantify  ;h“  expected 
background,  so  that  targat-background  contrasting  features  might  be  Identified, 
thus  Indicating  the  appropriate  specifications  for  detection  hardware. 

Field  measurements  of  radar  return  could  be  used  for  determining  natural 
background.  However,  this  can  be  quite  expensive  and  time-consuming.  Alterna- 
tively, one  can  use  computer  simulation,  which  can  be  much  cheaper,  recognizing 
that  there  Is  a possible  problem  of  model  accuracy. 

This  report  describes  the  development  of  a simulation  system  to  be  used  In 
the  calculation  of  radar  return  from  forests,  where  the  micro-physics  of  radar 
scattering  Is  left  open,  allowing  Improvements  In  accuracy  as  Improvements  are 
made  in  the  physical  model . 

1 2 3 

The  framework  of  the  system  Is  a computer  model  of  trees  and  forests  ’ * , 
which  Mathematical  Applications  Group,  Inc.  (MAGI)  has  adapted  for  use  In  the 
study  of  radar  signatures.  As  described  tn  the  sections  which  follow,  the  code 
accepts  as  Input,  the  type  of  tree(s)  desired,  (from  a current  library  of  four 
tree  types)  and  the  distribution  of  these  trees  within  a forest  area.  From 
these  data,  and  the  angle  the  Incident  radar  makes  with  the  downward  vertical 
direction,  the  code  determines,  as  a function  of  height  above  ground,  the 
angular  distributions  of  hits  with  leaf  and  twig  surfaces  and  the  angular  and 
size  distributions  of  hits  with  woody  branches  and  trunks.  The  essential  use 
to  which  these  results  can  be  put  In  the  future  Is  to  reduce  the  exceedingly 
complex  geometry  of  a forested  area  to  a set  of  distribution  functions  for  each 
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of  a number  of  horizontal  single-reflection  layers.  Subsequently,  these  dis- 
tributions will  permit  the  determination  of  reflection  and  transmission  matrices 
for  each  layer  - each  matrix  representing  scattering  from  a set  of  incident 
channels  to  a set  of  emergent  channels.  From  these  matrices,  a reflection 
matrix  for  the  entire  system  can  be  constructed. 

The  sections  which  follow  describe  the  new  code,  RADOT  (radar  on  trees), 
developed  by  MAG)  to  obtain  the  leaf,  twig  and  branch  distributions  referred 
to  above.  (The  second-stage  matrix  code  was  not  within  the  scope  of  the 
current  work.) 

RADOT  is  currently  operational  on  the  IBM  36O  and  CDC  6400  and  6600  sys- 
tems. The  CDC  6400  version  has  been  Installed  at  the  Engineer  Topographic 
Laboratory  (ETL). 
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2. 


SIMULATION  OF  A FORESTS!!  AREA 


2.1  Preliminary  Generation  of  the  Forest  Configuration 

Prior  to  the  RADOT  code  firing  and  tracking  simulated  radar  rays  through 
a forested  area  (In  order  to  obtain  angular  and  size  distributions) , It  must 
first  have  available,  as  Input,  a complete  description  of  the  forest  of  In- 
terest. This  consists  of  two  parts.  The  first  part,  described  In  Section  2.2, 
fixes  the  location  of  trees  within  the  forest.  In  addition  to  Its  location, 
each  tree  must  also  be  described  by  Its  tree  type  (e.g.,  pin  oak,  cottonwood) 
and  by  Its  size  and  orientation  relative  to  a pre-determlned  prototype  tree. 

The  second  part,  discussed  in  Section  2.3,  determines  the  description  of  each 
prototype  tree. 

2.2  Forest  Generator 

4 

Program  FORGEN,  developed  previously  by  MAGI  , utilizes  user  Input  and 
built-in  distribution  functions  to  place  "tree  boxes"  througnout  a forested 
area.  A tree  box  Is  an  Imaginary  rectangular  parallelepiped  (RPn) , perpendl- 

l 

cular  to  the  ground,  which  Just  fits  over  a tree.  The  tree  Itself  Is  one  of 
the  prototypes,  (see  Section  2.3).  scaled  and  reorientated  (In  the  horizontal 
plane)  by  the  FORGEN  code. 

FORGEN,  Itself,  Is  probably  a great  deal  more  flexible  than  required  by 
RADOT  users.  The  code  can  generate  tree  boxes  within  specified  areas,  (.long 
boundaries,  or  at  specific  Isolated  points.  Only  the  area  option  Is  needed 
for  RADOT.  Also,  FORGEN  can  construct  a forest  consisting  of  several  distinct 
"patches",  but  for  RADOT  one  patch  should  be  sufficient.  Nevertheless,  the 
complete  code,  along  with  Its  description,  has  been  delivered  to  ETL. 
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The  output  of  FORGEN  Is  a permanent  file  (or  magnetic  tape)  containing, 
for  each  tree  In  the  forested  area,  a description  of  Its  enveloping  tree 
box.  The  tree  box  data  are  grouped  first  by  forest  patch  (If  more  than  one 
patch  exists)  and  then  by  tree  type  (If  more  than  one  tree  type  exists). 

These  organized  data  are  referred  to  as  the  "forest  library". 

2.3  Tree  Generator 

5 

The  tree  generator  was  developed  previously  by  HAG  I In  order  to  repre- 
sent, with  a high  degree  of  faithfulness,  the  woody  structure  and  leaf  distri- 
bution of  a tree.  Based  upon  user-input,  the  tree  generator  explicitly  de- 
scribes, via  connected  truncated  right  angle  conas  (TRCs)  the  trunk  and  the 
primary  and  secondary  branches.  The  size,  shape,  location  and  orientation 
of  each  primary  or  secondary  branch  Is  Included  In  the  tree  description. 

The  actual  "construction"  of  a trae  Is  based  on  three  prototype  structures: 
the  primary  branch,  the  segment  of  the  primary  branch  and  the  secondary  branch. 

The  main  trunk  of  the  tree  Is  a stack  of  TRCs.  Copies  of  the  primary 
branch  (see  below)  are  hung  on  the  main  trunk,  as  directed  by  Input.  Copies 
of  secondary  branches  (also,  see  below)  are  also  permitted  on  the  mein  stem. 

The  primary  branch  consists  of  stacked  sections  of  primary  branch  seg- 
ments. Each  of  these  segments  consists  of  a TRC  on  which  copies  of  secondary 
branches  are  hung. 

The  prototype  secondary  branch  consists  of  a TRC  capped  by  twig  and  leaf 
structure.  This  structure  Is  resident  tn  the  model  (one  structure  per  tree 
type)  and  can  be  modified  tn  Its  three  dimensions  by  different  amounts. 
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Two  different  computer  models  for  each  tree  are  constructed  by  the  tree 
generator;  a high  resolution  model  and  a low  resolution  model.  These  differ 
only  in  the  trentment  of  twigs  and  leaves.  The  high  resolution  model  is  used 
only  for  representation  of  close-up  trees  for  HAGI's  earlier  code  which  gen- 
erates actual  pictures^.  (This  model  is  not  used  by  RADOT.)  In  the  high  re- 
solution model,  the  twig  ar.d  leaf  structure  Is  completely  deterministic.  The 
low  resolution  model  Is  used  for  representation  of  far-away  trees  for  MAG  1 1 s 
earlier  code  and  also  for  radar  tracking  In  RADOT.  In  the  low  resolution 
model,  the  twig  and  leaf  structure  Is  replaced  by  homogeneous  "clouds"  of 
twigs  and  leaves,  as  described  below. 

The  tree  generator  constructs  the  tree  box  about  the  tree  as  noted  In 
Section  2.2,  above,  This  box  Is  then  divided  tn  a regular  lattice  of  15x15x25 
cells;  the  division  by  "25"  being  in  the  vertical  direction,  The  tree  compon- 
ents are  then  sorted  Into  the  lattice  cells.  (Note:  A long  branch  can  span 
two  or  more  such  calls.)  For  the  low  resolution  model,  computer  experiments 
are  first  performed  with  the  high  density  model  to  obtain  the  proper  twig  and 
leaf  cloud  densities  for  each  cell,  These  densities  are  then  assigned  to 
each  cell.  Note  that  In  all  subsequent  ray  tracing  with  the  low  resolution 
trees,  twig  and  leaf  "hits"  are  determined  statistically,  based  upon  the  cloud 
densities  of  the  lattice  cells. 

To  date,  MAGI  has  processed  four  prototype  trees  in  both  the  high  and  low 
resolution  versions  . These  tree  data  have  been  forwarded  to  ETL.  The  pro- 
cessed trees  Include: 

1 . Virginia  Pine 

2.  Pin  Oak 

3.  Cottonwood 

k.  "Averaged"  Ft.  Belvoir  deciduous  (mostly  Gum)  tree, 
without  leaves. 

y— — 

The  prototype  trees  can  be  scaled  In  height  (and  other  dimensions)  by  proper 
Input  to  the  FORGEN  code. 
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The  collection  of  processed  trees  Is  referred  to  as  the  "tree  library". 


2.4  Organization  of  Forest  and  Tree  Data 

The  forest  and  tree  data  must  be  made  available  to  RADOT  on  logical 
unit  ID8  (see  Section  4.3).  The  order  Is  as  follows  (all  In  binary): 

1.  Forest  library  (output  of  FORGEN) 

2.  A single  number  flag  ■ 78130 

3.  The  tree  library  (output  of  the  tree  generator) 
containing  each  prototype  tree  referred  to  by 
the  forest  library. 

4.  A single  number  flag  ■ 15626 

The  apparently  peculiar  nature  of  the  two  flags  (second  and  fourth  Items) 
has  logical  meaning  only  for  MAGI's  earlier  code  but  was  carried  over  for  con- 
venience for  RADOT. 

It  Is  currently  up  to  the  user  of  RADOT  to  organize  the  forest  and  tree 
data  as  outlined  In  this  section. 


3.  operation  of  radot 


3 . 1 Problem  Specification 

After  the  forest  and  tree  data  are  organized,  as  outlined  In  Section  2.4, 
above,  the  user  supplies  some  additional  preliminary  Information  to  RADOT. 
These  Include: 

1.  Title 

2.  Card  punch  and  output  tape  options 

3.  Number  of  height  regions  In  which  data  are  to  be  recorded 

4.  The  elevation  above  ground  for  these  height  regions 

5.  The  number  of  angular  bins  (equally  spaced  tn  cosine  from 
+1.  to  -1.)  In  which  angular  distributions  are  to  he  accumu- 
lated. 

6.  The  number  of  diameter  bins  (equally  spaced  from  user-supplied 
maximum  diameter  to  zero)  tn  which  branch  distributions  are 

to  be  accumulated. 

7.  The  center  and  horizontal  dimensions  of  the  section  of  the 
forest  to  be  examined  and  the  angle  that  the  Incident  radar 
makes  with  the  downward  vertical. 

8.  The  nominal  number  of  radar  rays  to  be  "fired"  at  the  forest. 

This  Is  somewhat  altered  by  the  code,  as  explained  below. 

3.2  Preprocessing  by  RADOT 

RADOT  will  accept  the  specification  data  as  supplied  In  Section  3.1  and 
will  perform  the  following  preprocessor  operations:' 

1.  Construct  a focal  (firing)  plane  ebove  the  forest  section  as 
shown  in  Figure  1.  This  plane  will  be  oriented  so  that  radar 
fired  with  negative  Z and  positive  X direction  components,  passes 
perpendicularly  through  the  focal  plane. 
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downward  vertical 


direction  of  Incident  radar 


Foca 1 P 1 ane 


2. 


2 

Consider  an  N regular  lattice  of  points  on  the  plane,  so 
that : 

a.  N ■ 8J+1 

where  J Is  an  Integer  determined  by  condition  "b" 

b.  Hv/fi  where  M Is  the  nominal  number  of  rays  specified 
by  the  user  and  N Is  as  close  to  /FT  at  possible  while 
still  satisfying  condition  "b" 

2 

The  N lattice  points  are  uniformly  located  on  the  focal  plane,  and  are 
so  situated  that  radar  rays  passing  through  these  points,  will  Intercept  the 
entire  forest  area  (but  no  additional  area)  specified  by  the  user,  as  described 
In  Section  3.1 . 

A 

3.  RADOT  will  then  divide  the  focal  plane  Into  smaller  subsections 
(essentially  rectangles  of  8x8  points)  and  will  determine,  for  each  subsection, 
which  of  the  tree  boxes  supplied  by  FORGEN  could  be  Intercepted  by  radar  rays 
entering  that  subsection.  This  sorting  Information  Is  saved  and  used  subse* 
quently  by  the  ray  firing  routines  as  a computer  running  time  saving  device. 

3.3  Ray  Firing  and  Scoring 

RADOT  considers,  one  at  a time,  each  grid  point  In  the  focal  (firing) 
plane  and  "fires"  a ray  through  that  point  at  the  user-specified  angle  with 
the  downward  vertical. 

For  each  grid  point,  the  code  will  consider,  one  at  a time,  each  tree  box 
selected  as  a possible  Intercept  by  the  pre-sorting  routines  referred  to  In 
Section  3.2  The  ray  is  tracked  through  the  grid  point  (on  the  focal  plane) 
towards  a tree  box.  If  It  does  Intercept  a given  tree  box  further  examination 
and  possible  scoring  are  required,  as  outlined  In  the  following  paragraph. 
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If  It  misses,  the  code  examines  the  next  possible  tree  box.  This  process  then 
repeats  until  all  the  pre-sorted  tree  boxes  are  exhausted.  At  that  time,  the 
code  proceeds  to  Investigate  the  next  grid  point  * a process  which  continues 
until  all  grid  points  have  been  examined.  At  that  point,  control  ts  sent  to 
the  edit  routine,  as  described  In  Seetlon  3 • ^ » below. 

If  a ray,  fired  through  a particular  grid  point,  does  Intercept  a tree 
box,  then  further  examination  Is  required  to  determine  If  any  branches  or 
clouds  of  twigs  or  branches  will  Intercept  the  ray.  To  accomplish  this,  the 
code  considers  the  15x15x25  Internal  cell  array  of  the  tree  (see  Section  2.3). 
The  ray  Is  tracked  through  those  cells  which  It  encounters  along  Its  path 
(Inside  the  tree  box)  until  It  either  hits  the  ground  or  exits  laterally  from 
the  tree  box. 

For  each  cell  traversed,  the  coda  determines  If i , 

1.  a branch  ts  contained  In  the  cell,  and,  If  so,  If  the  branch 
Itself  Is  actually  hit. 

2.  the  twig  cloud  density  Is  greeter  than  zero. 

3.  the  leaf  cloud  density  Is  greater  than  zero. 

If  material  Is  encountered  due  to  any  one  of  the  Items  listed  above,  (or 
for  any  combinations  of  such  Items),  the  following  8-word  record  Is  placed  on 
an  "Interaction  file"  (logical  unit  IMT)  for  aach  of  the  encounters: 


1.  Type  of  encounter  (!.«.,  branch,  twtg  or  leaf-  I,  2,  or  3 respectively) 

2.  Score  (-1  for  branches;  - statistical 

estimate  for  twigs  or  leaves  * see  below) 

3-5.  Directional  Cosines  (of  major  TRC  axis  for  branches,  of  the 

perpendicular  to  the  surface  for  twigs  or 
leaves) 

6.  Height  of  Hit  (-actual  height  for  branches; 

■average  height  of  cell  box  for  twigs  or 

leaves) 


7. 

Diameter  of  Branch 
at  Point  of 
Encounter 

(applicable  for  branches  only) 

8. 

Ray  Number 

(not  currently  used  but  available  for 

future  statistical  analyses). 

For  twigs  or  leaves,  the  score  (Item  2)  Is  of  the  form 
I - exp (-Ed) 

where  d Is  the  distance  of  the  ray  through  the  cell  and  E Is  a macroscopic 
"cross-section"  (l.e.,  probability  of  a "hit"  per  unit  distance)  which  Is 
directly  related  to  the  twig  or  leaf  cloud  density  assigned  to  that  cell. 

After  writing  the  8-word  record  for  each  kind  of  encounter  In  the  cell,  the 
code  proceeds  to  the  next  cell  (If  any)  to  continue  the  tracking  process  as  If 
no  hit  had  been  made. 

Upon  the  completion  of  all  the  rays,  the  Interaction  tape  la  rewound  and 
sorted  In  the  edit  routine  (Section  lA), 


3 .*♦  Editing  of  Results 

Thfc  edit  routine  rewinds  the  Interaction  tape  (described  In  Section  3*3) 
and  reads  the  data  back,  event-by-event.  Each  event  Is  sorted  according  to 
Its  encounter  type  (l.e.,  branch,  twig  or  leaf)  and  Its  height  bln.  A counter 
for  the  appropriate  type-height  bln  Is  Incremented  by  the  event  score  (see 
Section  3*3).  In  addition,  a similar  Increment  Is  made  In  the  appropriate  type- 
height-angular  bln  - the  angular  bln  currently  beLng  described  only  by  the 
cosine  of  the  angle  between  the  upward  vertical  vector  and  the  score  direction 
cosines  (l.e.,  the  TRC  axis  for  branches;  the  normal  to  the  surface  for  twigs 
and  leaves).  The  angular  bln  structure  was  specified  by  the  user  upon  Input. 
Finally,  If  the  encounter  was  with  a branch,  an  additional  Increment  Is  made 
In  the  appropriate  helght-dlameter  bln  - the  diameter  bln  currently  being  de- 
scribed by  the  fraction  of  a maximum  user-speclf led  diameter. 

Upon  the  completion  of  all  sorting  and  Incrementing,  RADOT  will  then 
normalize  all  scores  to  a "per  ray"  basis  and  will  further  normalize  all  distri- 
butions so  that  their  Integrals  are  unity.  All  scores  are  then  printed  out  and 
options  are  available  to  punch  out  these  scores  or  to  write  them  on  output  mag- 
netic tape. 

Finally,  the  edit  routine  will  examine  the  total  scores  In  each  height  bln 
and  will  determine  and  print  out  the  subdivisions  of  each  bln  so  that,  assuming 
an  originally  uniform  height  distribution  of  scores  within  a given  height  bln, 
the  probability  of  a hit  In  any  given  subdivision  does  not  exceed  a fixed  prob- 
ability, Pj.  Currently,  (looking  ahead  to  a single  scattering  matrix  code), 

2 

P,  Is  set  to  0.1.  That  Is,  the  probability  of  two  or  more  hits  goes  as  MJ.l  or 
'v/.Ol  and  so  the  chance  of  a multiple  scatter  wlthl.n  a subdivided  height  bln  Is 
very  small . 
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4.  CODE  DESCRIPTION 


4.1  Routines 

The  following  paragraphs  provide  a brief  description  of  each  routine  of 
the  RADOT  code. 

MAIN  (RADOT) 

The  MAIN  program  Initializes  some  constants  and  reads  Input  data  describ- 
ing the  scoring  regions,  the  distributions  desired  and  the  types  of  trees  to 
be  Included. 

It  then  calls  CAMERA  to  locate  the  "firing  position"  and  the  focal  (l.e., 
"firing")  plane. 

Subsequently,  In  a loop  over  the  different  types  of  trees,  the  following 
routines  are  called: 

(1)  FORRET  - which  extracts  from  the  data  base  the  portlon(s)  of  the  pre- 
calculated forest  which  lie  In  the  specified  Investigation  area. 

■ 

(2)  PREPAR  - which,  In  order  to  speed  up  the  ray-t racing  In  the  latter 
stages  of  the  code,  examines  the  various  rectangular  areas  on  the  focal  plane, 
and  records  which  specific  trees  can  be  hit  by  rays  fired  from  the  firing  posi- 
tion through  that  rectangular  area. 

(3)  MAMAQE  - which  controls  the  ray-tracing  through  the  trees  and  the 
writing  of  Interaction  "hit"  data  on  the  Interaction  tape. 

When  all  tree  types  have  been  analyzed,  EDIT  Is  eat  led  to  edit  the  results. 
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CAMERA 


Subroutine  CAMERA  reeds  In  data  describing  the  forest  area  to  be  Investi- 
gated and  the  angle  of  the  Incident  rays  with  the  downward  vertical.  Based 
upon  this  Information,  CAMERA  computes: 

(1)  the  firing  position  for  all  rays, 

(2)  the  location  and  Inclination  of  the  focal  plane, 

, 4 

(3)  the  first  grid  point  (in  tha  focal  plane)  to  be  examined,  and 

(4)  data  which  will  subsequently  allow  subroutine  VPICTU  to  determine 
all  other  grid  points  In  the  focal  plane. 

EDIT 

Subroutine  EDIT  terminates  and  rewinds  the  interaction  tape  which  Is  de- 
signated as  logical  unit  "INTn.  EDIT  will  then  read  l?ack  each  Interaction 
event,  one  at  a time,  and  first  determine  the  height  region  of  the  event  and 
the  type  of  event  (l.e.,  solid  wood,  twigs  or  leaves)^  Tha  number  of  hits  by 
height  and  type  are  Incremented.  The  event  data  ere  then  further  evaluated  to 
determine  the  proper  polar  cosine  bln,  (of  the  axis  for  solid  wood,  and  of  the 
normal  vector  for  twigs  and  leaves),  and  the  appropriate  Incrementing  Is  per- 
formed. If  the  event  Is  with  the  solid  wood,  the  diameter  of  the  wood  Is  re- 
trieved from  the  data  and  the  appropriate  Incrementing  Is  performed. 

This  procedure  Is  continued  until  all  Interaction  events  have  been  re- 
trieved, examined  and  sorted.  At  this  point  all  results  are  printed  out,  In- 
cluding both  absolute  and  normalized  angular  and  diameter  distribution  date. 
Upon  option,  the  results  may  also  be  punched  out  or  written  on  an  output  file. 
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As  a final  step,  EDIT  will  examine  the  accumulated  data  In  each  large 
(input)  height  region  and  will  determine  the  layer  boundaries,  within  each 
height  region,  so  that  the  probability  of  a "hit",  per  Incident  ray,  Is  no 
greater  than  "Pj".  Note  that  P^  Is  set  by  an  executable  statement  In  the 
MAIN  code.  It  Is  currently  set  to  0.1  so  that  the  chances  of  a multiple- 
hit  goes  as  MO  which  Is  a small  enough  probability. 

FORRET 

FORRET  Is  called  by  the  MAIN  routine  once  for  each  type  of  tree  speci- 
fied by  the  user. 

FORRET  than  performs  the  following  operations,  utilizing  the  forest  data 
file  on  logical  unit  "IDB". 

(1)  Verifies  that  the  specified  tree  type  exists  on  the 
tree  library, 

(2)  Examines  each  patch  (sub-region)  of  the  forest  and 
determines  whether  that  patch  contains  trees  of  the 
type  currently  under  question.  If  It  does,  FORRET 
calls  INROUT  to  determine  whether  the  user-specified 
area  to  bo  examined  lias  In  the  patch.  If  It  does 
not,  the  routine  proceeds  to  the  next  patch.  If  It 

l 

does,  FORRET  determines  which  trees  (l.e.,  tree  boxes) 
lie  In  the  user-specified  area,  and  writes  the  tree 
box  Information  for  those  trees,  on  logical  unit  "IPD". 

Up  to  5000  tree  boxes  (of  a certain  type  of  tree)  may 
be  found  to  exist  within  the  user-specified  area. 
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ggbox 


GGBOX  Is  a routine  to  determine  If  • rey  originating  at  a specified 
point  and  aimed  at  a specified  direction  will  Intersect  a given  (tree)  box. 

If  nof  a logical  variable  "(HISS"  Is  set  to  "TRUE".  If  yes,  IMISS  Is  set 
to  "FALSE"  and  the  distance  to  the  box  Is  determined. 

GGTRC 

GGTRC  Is  similar  to  GGBOX,  except  that  It  Is  used  for  a truncated  right 
cone  (TRC) . 

INROUT 

INROUT  determines  If  any  portion  of  the  user-specified  aree  to  be  exam- 
ined lies  within  any  portion  of  a given  peteh  of  the  forest. 

LRAVT2 

LRAYT2  controls  the  tracking  and  scoring  (If  any)  of  a ray  through  a 
given  tree  box. 

Yha  tree  box  Is  divided  Into  a 15x15x25  cell  mesh  - with  the  25  sub- 
division being  along  the  height  vector.  Each  of  these  boxes,  according  to 
the  tree  type,  may  or  may  not  contain  solid  wood,  twigs,  or  leaves  - or  com- 
binations of  two  or  three  of  these  components.  The  solid  wood  Is  explicitly 
located  In  space;  whereas  twigs  and  leaves,  If  they  exist  In  a cell,  are 
assumed  to  be  uniformly  distributed  throughout  the  cell. 

The  ray  will  be  tracked  by  LRAYT2  from  the  point  where  It  enters  the 
tree  box  to  the  point  where  It  exits  - either  through  e side  or  by  striking 
the  ground. 

LRAYT2  calls  LRIGID,  which,  In  turn,  tracks  through  the  cells  until  one 
Is  found  that  ts  not  empty.  LRIGID  will  return  to  LRAYT2  the  number  of  solid 
wood  components  In  the  current  cell  (If  any)  and  the  Interaction  cross  sec- 
tion data  for  the  homogenized  twigs  (If  any)  and  the  homogenized  leaves  (If  any). 

' 

i 


16 


For  each  solid  wood  component,  LRAYT2  will  call  GGTRC  to  determine 
whether  an  actual  intercept  was  made.  If  so,  the  height  of  the  hit  (average 
height  of  the  path  through  the  TRC),  the  diameter,  and  the  normalized  direct- 
ion cosine*  of  the  TRC  axis  will  be  recorded,  along  with  other  Information, 
on  the  Interaction  tape  ("ItyT"). 

For  twigs  and/or  leaves  a score  of  the  type  [l-exp(-Ed)] , where  E Is  the 
macroscopic  cross  section  and  d Is  the  distance  through  the  cell,  will  be  com- 
puted for  E>0.  The  average  height  of  the  ray  through  the  cell  and  the  direct- 
ion cosines  of  the  normal  to  the  twig  or  leaf  will  also  be  computed.  These 
data,  among  others,  will  then  be  written  on  the  Interaction  tape. 

LRAYT2  saves  considerable  computer  time,  while  tracking  through  the  tree 
box,  In  the  following  manner.  If,  in  the  course  of  tracking  through  a specific 
cell,  GGTRC  determines  that  a certain  piece  of  solid  wood  (TRC)  Is  not  Inter- 
cepted by  the  ray,  a flag  Is  set  so  that  the  TRC  will  not  be  examined  again 
(for  that  ray)  for  any  other  call  In  which  It  may  appear. 

LR  IQ  I D 

LRIGID  tracks  from  cell  to  cell  within  a tree  box.  If  LRIGID  determines 
that  on*  or  more  TRCs  exist  In  a cell  and/or  homogenized  twig  or  leaf  material 
exists,  It  will  return  such  Information  back  to  Its  calling  routine,  LRAYT2. 

The  Information  about  what  Is  contained  tn  a given  cell  Is  retrieved  from  the 
CELL  and  CELL2  arrays  both  of  which  were  transferred  to  LRIGID  via  COMMON/TREE/. 
These  data  were  originally  set  by  the  portion  of  the  tree  cod*  (of  the  earlier 
work)  whleh  converts  the  tree  from  the  high- resolution  to  the  low-resol ut ton 
model . 

When  LRIGID  returns  to  LRAYT2,  a logical  variable  "EXIT"  Is  set.  If  the 
ray  Is  still  within  the  tree  box,  EXIT  ts  set  to  "FALSE".  If  the  ray  has  es- 
caped from  the  tree  box,  EXIT  Is  set  to  "TRUE". 
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MAMAGE 


MAMAGE  Is  called  by  MAIN,  for  each  tree  type,  to  control  the  tracking 
part  of  the  coda.  MAMAGE  reads  In  the  tree  data  from  logical  unit  "IPD" 
and  calls  VPICTU  to  begin  the  ray  tracing. 

PREPAR 

PREPAR  Is  called  by  MAIN,  for  each  tree  type.  Its  purpose  Is  to  pre- 
sort the  tree  boxes  (In  the  user-selected  part  of  the  forest)  so  that  only 
those  tree  boxes  appropriate  for  a given  rectangular  area  on  the  focal  plane 
(l.e.,  only  those  tree  boxos  which  can  be  Intercepted  by  a ray  from  the  firing 
point  aimed  through  the  given  grid  point)  need  be  considered. 

RANF 

Function  RANF  Is  a random  number  generator.  It  selects  a floating  point 
number  with  uniform  probability  between  0.  and  1. 

REED 

REED  Is  an  auxiliary  routine  called  by  PREPAR  to  read  an  array  from  a 
specified  logical  unit  In  the  minimal  amount  of  time. 

VPICTU 

VPICTU  selects,  one  at  a time,  the  grid  points  on  the  focal  plane  and 
starts  all  rays  towards  the  trees  at  the  user-specified  angle  to  the  vertical 
downward  direction.  It  then  calls  VRAYTR  to  continue  tho  tracking  process. 
VRAYTR 

VRAYTR  considers  the  focal  plane  grid  point  In  question  and  retrieves, 
from  logical  units  3 and  10,  the  data  for  the  appropriate  tree  boxes  (l.e., 
those  selected  earlier  by  PREPAR).  It  then  calls  VRAYT1  to  continue  the 
tracking  process. 
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VRAYT1 


VRAYT1  considers  one  at  a time,  the  tree  boxes  (for  a given  focal  plane 
grid  point)  as  retrieved  by  VRAYTR.  For  each,  VRAYT1  calls  GGBOX  to  deter- 
mine tf  the  ray  does  actually  Intercept  the  tree  box.  If  not.  It  proceeds 
to  the  next  tree  box.  If  It  does,  control  Is  turned  over  to  LRAYT2  for  the 
actual  tracking  through  the  tree  box.  When  all' possible  tree  boxes  have  been 
examined,  VRAYT1  returns  control  to  VRAYTR. 

WRIT 

WRIT  ts  an  auxiliary  routine  called  by  PREPAR  to  write  an  array  on  a 
specified  logical  unit  In  the  minimal  amount  of  time. 

4.2  Code  Structure 

Figure  2 displays  the  organization  of  the  RADOT  code  and  the  labeled 
commons  utilized  by  each  routine. 

4.3  Logical  Unit  Assignments  ; 

The  card  reader  and  printer  are  assigned  units  IN  and  OUT,  respectively. 
The  forest  and  tree  data  must  be  supplied  from  unit  IDS  and  the  Interaction 
tape  Is  unit  INT.  An  Important  Intermediate  "scratch"  file  Is  unit  I PD* . 
Currently,  these  units  are  assigned  by  executable  statements  In  the  MAIN 

i 

routine  as  follows: 

1 08*1 
IN-5 
OUT-b 
I NT-1 4 
IPO-15 

In  addition,  units  3,  10,  11,  12  and  13  are  use$  as  Intermediate  scratch 
files.  Separate  unit  assignments  can  ba  made  by  the  user  (or  the  optional 

it 

punch  and  output  flies  . 

NOTE:  If  desired,  file  IPD  can  be  rewound  by  the  edit  routine  and  used  for 
tape  (or  disk)  output.  Tape  7 Is  usually  used  for  punching  cards, 
tf  that  option  Is  desired. 

I ' 
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FIGURE  2 - Cod*  Structure  of  RADOT,  with  1ab*1*d  COMMONS 
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4,4  Error  STOP  Messages 

Various  error  conditions  will  cause  the  RADOT  code  to  terminate  with 
numbered  STOP  statements.  In  such  cases,  the  current  routine,  itself,  pro- 
vides information  as  to  the  cause  of  the  STOr1  condition.  As  a user  aid,  the 
locations  of  the  various  STOP  conditions  are  now  provided: 


ROUTINE 

STOP  Number (s) 

MAIN  (RADOT) 

11 

CAMERA 

101,  102 

EDIT 

106,  107,  113 

F0RRET 

20,  21,  22 

LRAYT2 

40,  42,  44 

PREPAR 

60 

4.5  Input  Description 

The  following  Is  a complete  description  of  the  input  required  by  RADOT: 

1.  TITLE  (Format  20A4) 

2.  NHAC, I PUN, I TAP  (Format  3110) 

NMAC  ■ Number  of  height  layers  for  which  distributions 
shall  be  determined.  (<_1 00) 

IPUN  ■ Punch  option;  -0  for  no  punch,  ■ punch  unit  assignment  for  punch. 
ITAP  - Tape  output  option;  «0  for  no  output  tape,  - tape  unit  assignment 
for  output  tape. 

3.  MITE  (Format  6E12.4) 

NMAC+i  boundary  values  (feet)  of  the  leyers. 

HITE(l)  must  not  be  lower  than  the  top  of  any  tree. 

KITE (NMAC+l)  must  be  0. 
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4.  NTH  ET , MO  I AH , D I AMAX 


(Format  2110, E12. 4) 


NTHET  ■ number  of  equally  spaced  polar  angle  bins  for 
distributions* 

NDIAM  ■ number  of  equally  spaced  diameter  bins  for 
distributions* 

D I AMAX  ■ maximum  diameter  of  a tree  (feet) 

5.  XO,YO,XL,YL,COSTH  (Format  5EI2.4) 

XO  ■ center  of  forest  In  X-dlrectton  (feet) 

'<  Y0  ■ center  of  forest  In  Y-dlrectlon  (feet) 

XL  ■ X-length  of  forest  (feet) 

YL  - Y-length  of  forest  (feet) 

COSTH  ■ polar  cosine  of  rays  from  vertical  downward  direction 

6.  IFR  (Format  110) 

lilt 

IFR  ■ spectflc  forest  patch  to  be  examined 
■ 0 neglect  this  option 

7.  MRAYS  (Format  110) 

MRAYS  - nominal  number  of  rays  to  be  fired 

8.  JFIRST, JLAST  (Format  2110) 

JFIRST  ■ tree  no.  for  first  tree 
JLAST  ■ tree  no.  for  last  tree 


NMAC*NTHET<2000}  NMAC*NDIAM<2000 

If  the  forest  contains  more  than  one  pateh  and  consideration  of  only 
one  patch  ts  desired. 


JUJUX  [ 

Use  1,2, 3,4  for  Virginia  Pine,  Pin  Oak,  Cottonwood,  Ft.  Belvolr- 

averaged  Deciduous  (without  leaves),  respectively.  j 

I 

I 

I 

i I 


4.6  Output  Description 

4.6.1  Organization  of  Output 

The  printed  output  consists  of  three  parts: 

1 . Print-beck  of  Input 

2.  Intermediate  results 

3.  Results  of  ray  firing 

These  three  parts  are  discussed  In  the  saettons  which  follow. 

4.6.2  Print-Back  of  Input 

RADOT  prints  back  for  the  user  the  Input  as  specified  In  Section  4.5, 
above.  In  addition  to  the  nominal  number  of  rays,  the  actual  number  of 
rays  Is  printed  out,  as  wall  (see  Section  3 .2). 

4.6.3  Intermediate  Results 

Some  of  the  Intermediate  results  derived  by  RADOT  prior  to  ray  firing 
are  printed  out  as  extra  Information  for  the  user.  These  Include  the 
following: 

1.  Firing  Position  - the  nominal  location  of  a firing  position 
used  to  pre-sort  the  tree  boxes. 

2.  Start  of  Focal  Plane  - the  first  point  on  the  focal 
(firing)  plane. 

3.  Increments  - As,  and  Ay;  used  In  finding  other  points 
on  the  focal  plane. 

4.  Forest  Description  - Including  for  each  patch  In  the 
forest,  (there  may  be  more  than  one): 

a.  pateh  number 

b.  user-selected  patch  (If  any) 

c.  number  of  trees  by  tree  type. 
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Under  Item  "c",  the  tree  type  Is  a two  digit  number,  where  the  first  digit 
Is  the  tree  type  end  the  second  digit  Is  I or  2,  depending  upon  whether  the 
forest  was  described  by  the  high  or  low  resolution  model,  respectively.  Al- 
though RADOT  uses  the  low  resolution  tree  model  only,  It  can  pick  up  and 
utilize  a forest  description  (from  FORGEN)  with  the  high  resolution  tree  box 
since  the  tree  boxes  for  both  resolution  models  are  Identical. 

5.  Number  of  Trees  before  Sorting  - the  number  of  trees 
selected  as  potentially  being  hit  by  rays  fired  through 
the  selected  points  on  the  focal  plane. 

6.  Tree  Sorting  - the  sorted  end  selected  tree  data  ere  stored 
on  scratch  files  II,  12  and  13.  This  edit  Indicates  how 
many  trees  are  on  each  file.  Since  a tree  can  be  entered 
on  more  than  one  file,  the  sum  of  the  trees  listed  here 
can  exceed  the  number  given  In  Item  5. 

7.  Tree  Assignments  (repeated  for  flies  11,  12  and  13)  - 
the  number  of  trees  for  each  fopal  plane  area  asstgned  to 
this  file  Is  given,  as  well  as  the  tree^numbers  In  the 
table  which  follows  as  Item  8. 

8.  T-ae  Origins  (repeated  for  files  11,  12  and  13)  " 

the  tree  number  and  the  tree  box  x and  y origin  for  each 
tree  assigned  to  this  file. 

9.  Number  of  Trees  after  Sorting  - some  trees  originally 
selected  In  Item  5 may  have  been  eliminated  from  further 
consideration.  Hence  Item  9 is  not  necessarily  Identical  to 
Item  5. 

10.  Forest  Complete  through  Tree  Type  - printed  out  for  each 
completed  tree  type. 
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Results  of  Ray  Firing 

The  final  edit  is  printed  out  for  the  user,  as  described  In  Section  3.4. 
Punched  out  results  and  results  on  an  output  disk  file  (or  magnetic  tape)  are 
available  as  well  as  options. 

The  sample  problem  shown  In  Section  5,  Illustrates  the  format  of  the 


output. 


5.  TEST  PROBLEM 


5. 1 Description  of  Problem 

The  demonstration  problem,  described  In  this  section,  consists  of  rays 
fired  at  a patch  of  a forest  previously  generated  at  MAS  I by  the  F0R6EN  code. 

This  patch,  (number  1 of  the  entire  forest),  consists  of  tree  types  I 
and  2 In  the  tree  library,  which  are  the  Virginia  Pine  and  the  Pin  Oak,  re- 
spectively. As  In  all  runs  with  RAOOT,  the  low-resolution  tree  models  ere 
used. 

The  trees  are  randomly  mixed  In  the  forest  with  the  pine  tree  boxes  and 
oak  tree  boxes  covering  501  and  40*  of  the  forest  surface,  respectively.  The 
pine  trees  are  5. 0+0. 5 feet  In  height  while  the  oak  trees  are  12.0+0.2  feet 

The  number  of  rays  that  are  tracked  Is  1089  (In  a 33x93  grid).  They  are 
fired  at  a segment  of  the  selected  forest  patch  which  extends  In  the  "x-dlrect- 
ton"  from  360  to  420  feet  and  In  the  "y-d fraction",  from  700  to  760  feet.  There 
are  86  pine  trees  and  57  oak  trees  in  this  rectangular  segment.  The  rays 
Intercept  the  forest  at  an  angle  of  25.842*  (t.e.,  cosine  of  0.900)  with  the 
downward  vertical. 

Scores  are  accumulated  In  3 altitude  bins  with  height  boundaries  of  12.3* 
8.2,  7.2  end  0.0  feet. 

: 

Six  polar  bins  are  used  for  the  cosine  distributions  of  the  normal  vectors 
to  the  surfaces  of  the  twigs  and  leaves  and  for  major  axes  of  the  woody  seg- 
ments. Five  diameter  bins  are  used  for  the  woody  segments. 

The  problem  ran  for  84  seconds  on  an  IBM  360/65  computer. 

5.2  Input 

The  page  which  follows  displays  the  Input  to  the  test  problem. 
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FIGURE  3 - Input  to  Test  Problea 


5.3  Results  of  Ray  Firing 

The  pages  which  follow  display  the  output  edit.  (This  format  may  ba 
altered  by  ETL  In  the  future.) 

Note  that  for  the  distributions,  the  first  sat  of  data  are  the  absolute 
results  (hlts/ray)  and  the  second  set  of  data  are  the  results  normalized  to 
unity. 

Also  note  that  the  total  number  of  hits  between  heights  of  7.2  faet  and 
0.0  feat  Is  0.2879  par  ray.  Since  this  exceeds  the  maximum  allowed  value  of 
0.1,  RADOT  breaks  this  range  Into  Intermediate  ranges  with  the  Intermediate 
boundaries  being  printed  out.  For  the  other  two  height  layers  the  number  of 
hits  per  ray  are  less  than  0.1  and  so  no  further  subdivision  Is  required. 
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FIGURE  ^ - Output  of  Test  Proble* 
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